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Abstract 

We explore the energy spectrum of a non-relativistic particle bound in a linear finite 
range, attractive potential, envisaged as a quark- confining potential. The intricate tran- 
scendental eigenvalue equation is solved numerically to obtain the explicit eigen-energies. 
The linear potential, which resembles the triangular well, has potential significance in 
particle physics and exciting applications in electronics. 
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1 Introduction 

A challenging problem in particle physics in recent years is that of quark confinement. 
It is presently known that mesons are not elementary particles, but are composed of 
quarks, as are the nucleons. 

In literature, several approximation methods are available relating to quark confine- 
ment. The lattice model [1] suggests that at large distance between quarks, the interac- 
tion increases approximately linearly with separation. The bag model, where quarks and 
gluons are confined in a bag, is not suitable for calculating the hadronic properties of 
heavy quarks or in computing the energy levels of excited states. String model, on the 
other hand, proposes quark- ant iquark pair at the ends of an open string and creation 
of quark-antiquark pair when the string breaks. In recent years, potential models [2] are 
best justified theoretically to describe heavy quarkonia and seem to be most powerful in 
calculating the static properties. 

Several authors [3 - 6] have addressed the bound states of various kinds of linear 
potential. Chiu [7] has examined the quarkonium systems with the regulated linear 
plus Coulomb potential in momentum space. Deloff [8] has used a semi-spectral Cheby- 
shev method for numerically solving integral equations and has applied the same to the 
quarkonium bound state problem in momentum space. 

Rao and Kagali [9 - 11] have investigated extensively the bound states of both spin 
and spinless particles in a screened Coulomb potential, having linear behaviour near the 
origin. In the present work, we propose a finite, short-ranged linearly rising potential, 
envisaged as a quark-confining potential and explore the non-relativistic bound states. 

2 The Schrodinger equation with the linear potential 

Several attempts have been made to study the meson spectra using the non-relativistic 
Schrodinger equation with a linear potential. Intuitively, we construct a simple linear 
rising, finite range potential of the form [12] 



in which the well depth V and range 2a are positive and adjustable parameters. The 
linear potential with its boundary regions is illustrated in Fig.l and owing to its shape, 
this potential could also be called the trianglular potential well. 

Obviously in regions I and IV, the particle is free and the allowed solutions of the 
free particle Schrodinger equation are 




(1) 



a 




(2) 
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ij}\{x) = C 6 e ax a < x < oo, 
consistent with the requirement if)(x) vanishes as \x\ — > oo. 

Here a 2 - 



(3) 



- 2 i s implied. Since -E < for bound states, a is positive. To discuss 

h 

the nature of the soution within the potential region, —a < x < a, we insert the potential 
described in Eqn.(l) in the celebrated Schrodinger equation and obtain 



d 2 ip 2m 
dx^ + l T 



x\ 



E + Vq[1 — — 



ip(x) = 0. 



(4) 



x 



E - V 

Setting ~= U and defining E ~ ^2^2 ma 2 an< ^ ^° = /^/2ma 2 W6 °^^ a ^ n ^ ne dimensionless 
form of the Schrodinger equation 



<*V(y) 

dy 2 



+ 



E + V (l-y) ^{y) = 0, 



which may further be written as 

d 2 ^ 
dy 2 



Aytp + Bi\) = 0. 



(5) 



(6) 



The constants A — V and B — E + V are also dimensionless. Introducing an auxiliary 
function 

(7) 



w = (y - — ^ 



yeilds 



d 2 ^> 



— = 0. 



(8) 



The solutions of this differential equation are the well-known Airy functions Ai(w) 
and Bi(w) [13], having oscillatory and damping nature. 

The Eigenvalue Equation 

The admissible solutions in the four regions, consistent with physical reality, are 



— oo < x < —a 



(9) 



il> 2 (x) = C 2 Ai{w) + C 3 Bi(w) 



a < x < 



(10) 



^ 3 (x) = C A Ai(-w) + C b Bi(~w) 



0< x<a 



(11) 



i) A (x) = C 6 e~ 



a < x < oo 



(12) 



where C\ to C 6 are the normalisation constants. Imposing on the solutions in equations 

(9) to (12) the requirements that ip and ^ be continuous at the origin and also at the 
potential boundaries (x = ±a) leads to the eigenvalue eqnation. 



At x = -a, ^i(x) = i/j 2 (x) and -jj^=-^ 



This leads to 



C ie - aa = C 2 Ai( Wl ) + C 3 Bi(w) 



Ai 

aC ie - aa = — [C 2 Ai'( Wl ) + C 3 Bi'( Wl )\ 
a 



where 



Wl = As (-1 - 



B 



On simplification, we obtain 



a = 



P Ai'iwx) + Bi'(w 1 ) 



P Ai(wi) +Bi(wi) 



C 2 



where P = ~q^ ■ Similarly, the continuity condition at x = demands 

*(»)=*(*) and*=* 
from which we obtain 

C 2 Ai(w ) + C 3 Bi(w ) = C A Ai(-w ) + C 5 Bi(-w ) 
C 2 Ai'(w ) + C 3 Bi'(w ) = C 4 Ai'(-w Q ) + C 5 Bi'(-w ) 



with 



As before, 



w = As 



-B 
~A 



P Ai'(w ) + Bi'jwo) _ Q Ai'(-wo) + Bi'(-w ) 
P Ai(w ) + Bi(w ) ~ Q Ai(-wo) + Bi(-w ) 



where Q = qT is another constant. Adapting similar procedure at the boundary x 
demanding ip 3 (x) = i^a{x) and ^T = ^T one would on similar grounds obtain 



—a 



Q Ai'(-w 2 ) + Bi'(-w 2 ) 
Q Ai(-w 2 ) + Bi(-w 2 ) 



with 
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It is worthwhile mentioning that the arguments of the Airy function wq, Wi, and w 2 
are dependent both on the energy as well as the potential and are related by the simple 
equation 



w = 



w 1 + w 2 



(23) 



It is straightforward to check that 

_ Bi( Wl ) - Ak Bi'jwx) 
As Ai^wx) - (3 A^wx) ' 



(24) 



and 



Q = - 



P Bi{-w 2 ) + As Bi'(-w 2 ) 



(25) 



P Ai(-w 2 ) + As Ai'(-w 2 
where P = aa is implied. 

Formally on eliminating P and Q in Eqn.(20) we obtain the eigenvalue equation as 

'{pBi (oji) - AsBi' (wi)} Ai' (uq) + {AsAi' M - pAi ( Wl )} Bi' (uj ) 
[pBi (ujx) - A*Bi' (u x )}Ai (cu ) + {A*Ai' (wi) - PAi M} Bi (cu ) 



(26) 



{pBi(-uj 2 ) + AlBi'{-uj 2 ) 


}Ai'(-cj )- 


[A\Ai'{ 


-u 2 )+pAi(-u 2 ) 


}Bi'(-u ) 


[pBi(-uj 2 ) + A\Bi' (-u 2 ) 


}Ai (-w ) - 


{A^Ai' 


'-u 2 )+pAi(-u 2 ) 


}Bi (-u ) 



This intricate and fairly complicated transcendental eigenvalue equation involving 
the Airy function and its derivatives is solved both graphically and numerically using 
Mathematica[14] . The real roots, which correspond to the eigenenergies, are listed in 
Table 1 for a typical value of the range parameter (a). Energy (E) and well-depth (V ) 

are both expresed in units of 2 ma 2 ■ 



3 Results and Discussion 



One of the distinctive characteristics of quantum mechanics, in contrast to classical 
mechanics, is the existence of bound states corresponding to discrete energy levels. It 
is well-known in quantum mechanics that bound states exist for all atractive potentials, 
the exact number depending on the specific form of the potential and the dimensionality 
of the space. 

More specifically, as is seen from the spectrum of energies listed in Table 1, for a 
finite range of the potential, deeper wells admit excited state energies, consistent with the 
wisdom of quantum mechanics. Such studies, apart from being pedagogical in nature, 
are potentially exciting and significant as it is concerned with quark confinement. 
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Quantum chromodynamics, which governs the quark-antiquark interaction is widely 
accepted as a good theory of strongly interacting particles. One can explore the hadronic 
properties by investigating the bound states of quarks. Our investigation concerning the 
linear potential is seeming interesting and can be regarded as a model to describe the 
quarkonia. 

Further, the linear potential well or in other words, the triangular well has potential 
applications in electronics. Interestingly, in many semiconductor devices, it is believed 
that electrons are confined in almost triangular quantum wells [15]. Examples of such 
devices are Si MOSFETs (Metal Oxide Semiconductor Field Effect Transistors) widely 
used in digital applications and GaAs/AlGaAs MODFETs (Modulation Doped Field 
Effect Transistors) used for high speed applications. The bound states of the linear 
potential is a subject of renewed interest and intensive research and we have extended 
the study of this naive potential to the relativistic domain, which will be reported shortly. 
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Table 1 



Eigenenergies of a non-relativistic particle in 

linear potential 

(a = 1A) 
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0.50 


- 0.1122377 






0.60 


- 0.1451770 






0.70 


- 0.1790694 






0.80 


- 0.2134967 






0.90 


- 0.2480162 






1.0 


- 0.2828516 






2.0 


- 0.6121732 






3.0 


- 0.9072505 






4.0 


• 1.2056918 






4.28 


- 1.2984370 


- 0.0000991 




5.0 


- 1.5770587 


- 0.1653382 




10.0 


- 5.8335771 


- 1.3641773 




20.0 


-16.4420518 


- 2.3126501 




20.62 


-17.0607820 


- 2.3662130 


-0.0000343 


25.0 


-21.2732190 


- 3.6112921 


-2.4865189 


30.0 


-25.7497380 


- 9.0650370 


-2.9501500 


35.0 


-29.8836000 


-15.0604510 


-3.2815680 


40.0 


-33.7181130 


-21.1744464 


-3.5904940 



